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ABSTRACT. We present a class of graviton-dilaton models in 
which a homogeneous isotropic universe, such as our observed 
one, evolves with no singularity at any time. Such models may 
stand on their own as interesting models for singularity free cos- 
mology, and may be studied further accordingly. They may also 
arise from string theory. We discuss critically a few such possi- 
bilities. 
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1. General Relativity (GR) is a beautiful theory and, more importantly, 
has been consistently successful in describing the observed universe. But, 
GR leads inevitably to singularities. In particular, our universe according 
to GR starts from a big bang singularity. We do not, however, understand 
the physics near the singularity nor know its resolution. The inevitability 
and the lack of understanding of the singularities point to a lacuna in our 
fundamental understanding of gravity itself. By the same token, a successful 
resolution of the singularities is likely to provide deep insights into gravity. 

In order to resolve the singularities, it is necessary to go beyond GR, 
perhaps to a quantum theory of gravity. A leading candidate for such a 
theory is string theory. However, despite enormous progress string theory 
has not yet resolved the singularities 0. 

Perhaps, one may have to await further progress in string theory, or a 
quantum theory of gravity, to be able to resolve singularities. However, 
there is still an avenue left open - resolve the singularities, if possible, within 
the context of a generalised Brans-Dicke model, and then enquire whether 
such a model can arise from string theory, or a quantum theory of gravity. 

Brans-Dicke theory is special for more than one reason. It appears nat- 
urally in supergravity and Kaluza-Klein theories, and also in all the known 
effective string actions. It is also perhaps the most natural extension of GR 
[0], which may explain its ubiquitous appearance in fundamental theories: 
Dicke had discovered long ago the lack of observational evidence for the prin- 
ciple of strong equivalence in GR. Together with Brans, he then constructed 
the Brans-Dicke theory ^ which obeys all the principles of GR except that 
of strong equivalence. 

Brans-Dicke theory contains a graviton, a scalar (dilaton) coupled non 
minimally to the graviton, and a constant parameter uj. GR is obtained 
when u) = oo. In the generalised Brans-Dicke theory, commonly referred to 
as graviton-dilaton theory, u> is an arbitrary function of the dilaton [0, 
The graviton-dilaton theory contains an infinite number of models, one for 
every function very likely including all the effective actions that may arise 
from string theory, or a quantum theory of gravity. 

This perspective opens up an avenue in resolving the singularities: to 
find, if possible, the class of graviton-dilaton models which are singularity 
free. One may then study whether they can arise from string theory, or 
a quantum theory of gravity. On the other hand, such models may stand 
on their own as interesting models for singularity free cosmology and one 
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may then study their imphcations in other cosmological and astrophysical 
contexts, which are hkely to be novel 0. Such studies are fruitful and are 
likely to provide valuable insights. 

Accordingly, in this letter, we present a class of graviton-dilaton models 
where the function uj satisfies certain constraints. We show that a homoge- 
neous isotropic universe, such as our observed one, evolves in these models 
with no big bang or any other singularity. The time continues indefinitely 
into the past and the future, without encountering any singularity. These 
models also satisfy the observational constraints imposed by solar system 
experiments. We then discuss critically a few possibilities of such a model 
arising from string theory. 

2. Consider the following graviton-dilaton action 



S = J d'x^g i-xR + ^ (Vx)') + Sm{M, g,.) , (1) 



where is Newton's constant, x ^ is the dilaton, lj{x) is an arbitrary 
function characterising the theory, and Sm is the action for matter fields, 
denoted collectively by M. and coupling minimally to g^^. The "matter" is 
assumed to be a perfect fluid with density p and pressure p, related by p = 7p 
where — 1 < 7 < 1. The value of 7 indicates the nature of the "matter". For 
example, 7 = 0, |, or 1 indicates that the "matter" is dust, radiation, or 
massless scalar field respectively. If x has a finite range < x < Xmax < 00, 
as is the case in the present model, then the factor Xmax can be absorbed 
into On and the range of x can be set to be < x < 1. Also, we will work 
in the units where Gn = 

In the model we propose, derived originally in a different context in [^], 
X has a finite range, set to be < x < 1, and the function fi(x) = 2co'(x) + 3 
is required to satisfy the following constraints only [0: 

(i) nix) > 0, < X < 1 (ii) ^^(0) = ^0 < ^ 

(iii) limfi = fii(l-x)-'", ^<«<1 

(iv) -j-^ = finite V n > 1, < x < 1 ■ (2) 

Qq > and r2i > above are constants. For the purpose of illustration 
here, we further take ^{x) to be a strictly increasing function. However, 
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this constraint is not necessary. In fact, as will be shown elsewhere, the 
constraints (§) can also be relaxed further without affecting the results. 

Consider a homogeneous isotropic universe. Then, x = xi^) ^i^nd the line 
element is given by ds^ = —dt^ + e^"^*-*-* (dr^ + r'^dS'^)-, where ^5*2 is the 
line element on an unit sphere, and the curvature index k is set to zero. 
(The results, however, are valid for k = ±1 also.) The equations of motion 
following from the action (|l|) reduce to (see [Q) 

= + (3) 




2X 

P = Poe-'"*"'^ (5) 

where p = 7p is used, upper dots denote ^-derivatives, and po > is a 
constant. The square roots are to be taken always with a positive sign. 
Integrating equation (H) once, we get 

m = ^ {cr{t) + c) , a{t) ^ ^ / dt— ^ (6) 



and c is a constant. Dividing (|) by - and substituting (||) for x, we get 

Equations (|]) - (^ will prove useful in reading off the essential features of 
the evolution. 

3. We need to determine whether the evolution of the universe in the 
present model is singular or not. Directly obtaining explicit solutions to 
(0)"®) with Vt arbitrary, is not possible under general conditions pertaining 
to the observed universe. (For special cases, see ||^.) Nevertheless, by a 
careful analysis, it turns out to be possible to read off all the features of the 
general solutions. We will show this below. Particularly, we will show that 
a homogeneous isotropic universe, such as our observed one, evolves in the 
present model with no big bang or any other singularity. 

First, an observation. For the singularities to be absent, all curvature 
invariants must be finite. A sufficient condition for this, proved in the Ap- 
pendix, is that the quantities in (^) be all finite. However, some or all of 



4 



these quantities can potentially diverge when and/or x vanish or reach a 
critical value, or when Q ^ oo. Hence, we pay particular attention to these 
special points. 

To proceed, we need initial conditions. We start with an initial time 
ti, corresponding to a temperature > 10^^ GeV such that GUT symmetry 
breaking, inflation, and other (matter) model dependent phenomena may 
occur for t > only. The initial conditions at tj, relevant to our universe, 
are A{ti) > and xi^i) > which correspond to an expanding universe and 
a decreasing effective Newton's constant (= ig~)- The choice of ti means 
that, for t < ti, the universe is dominated by "matter" with 7 > | for as 
long as the scale factor is decreasing. We will first analyse the evolution 
for t < ti, which is when the big bang singularity arises in GR, and then for 
t>ti. 

3a. t < ti : Now, A{ti) and x{U) and, hence, the constant c in (j^) and 
^(tj) are positive. Therefore, K{ti) > 1 in (|^). As t decreases below ti, the 

quantities in K vary as follows: A decreases, since A > 0. The integral in (|^) 
is negative since t < ti, its prefactor is < 0, since A decreases and therefore 
7 > |. Therefore ait) is non decreasing. Equation (|^) then implies, since A 
decreases, that > x{ti) > 0. Hence, x decreases. In turn, fl decreases as 
X decreases, which follows from the properties of Q. The net consequence, 
therefore, is that K decreases monotonically as t decreases below tj. 

Note that the lowest value of X is ^ < |, achievable at x = 0. Since 
K{ti) > 1 initially, it follows from the above discussion that there exists a 
time tm < ti where K{tm) = 1 with x{tm) > 0. Then, ^{tm) = implying, 
since x(tm) > 0, that A{tm) = 0. Furthermore, since K{t) remains finite and 
x{t) remains non zero for tm ^ t < ti, it follows from (|^) that A{tm) = finite, 
implying that the scale factor e^'-*™-* is non zero. Thus, for t < ti, the universe 
continues to shrink and reaches a minimum non zero size at t = tm- (That 
this critical point of v4 is a minimum is physically obvious; it can also be seen 
from the positivity of A{tm)-) The precise values of tm and of A, x, and Q at 
tm are model dependent and nothing further can be said about them, except 
that fl{tm) < 3 (equality iff po = 0) since K{tm) = 1- 

However, the above information suffices for our purposes. For tm t < ti, 
X and remain non zero. Hence, it follows from (|^) that x{t) is finite. There- 
fore, the quantities in (|12|) all remain finite, implying that all the curvature 



invariants are also finite. Thus, there is no singularity for tm ^ t < ti |T0 
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3b. t < tra '• As t decreases below tmiA. becomes negative and, hence, 
A increases. The integral in is negative as before. However, its prefactor 
can be positive, since A increases now and therefore 7 can be < |. Hence, 
a{t) can be negative for t < tm- 

First consider the alternative that a{t) + c > for all t < t^. Then x > 0, 
implying that x continues to decrease for all t < t^n- Thus, singularities may 
arise as x ^ 0. In this limit, ^ il^. In obviously, 
either <^ 1, or 0{1), or ^ 1. Considering each of these possibilities in turn, 
and solving consistently equations (0) first and then (H), yields the unique 
solution as X — 0: 

forl^oT^^: e^-^^=(- + b)\ X = Xo (- + bV (8) 
6 Km / \m / 

for ^0 = ^ : e^-^« = e-f , X = Xoe"* , (9) 

where Aq, xo > 0, a > 0, and b are (model dependent) constants, and 

3 — vSHo —2 . . 

n = — , , m = ; . (10) 

3(1 - Vsn^) 1 - yMh 

Consistency of these solutions requires that 7 > (see below for 7 < 0) and 
\/3ilo > ^j3^. Note that f2o = | is perhaps the best choice, which is also 
significant from another point of view, as discussed at the end. 

Let f2o > I. Then, m > in (^0|). Since x ^ 0, it follows from (||) that 
t — > — ^ = finite, implying that x vanishes at a finite time in the past. Also, 
the quantities in (|12|), for example ^, diverge, implying that the curvature 
invariants, including the Ricci scalar, diverge. Thus, for fig > |) there is a 
singularity at a finite time in the past. 

Let fio < |, which is the case in our model, see (H). Then, m < in ([T0|) . 
Since x 0, it follows from (|^) and (|^) that t — > —00. Evaluating now the 
quantities in (p!2D, it can be seen that they are all finite, implying that all 
the curvature invariants are finite. Thus, for fio < there is no singularity 
at any time in the past. 

Now consider the remaining alternative that a{t) + c is not positive for 
all t < tm which will be the case, for example, for 7 < 0. Then, a{t) + c and, 
hence, x must vanish at some time t = tr, with xi^r) non vanishing. (The 
case of vanishing xi^r) is similar to the one discussed above.) Clearly, as 
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follows from equations (H), and (|^, all the quantities in (|T2|) remain 

finite at tr and, hence, there is no singularity at t = tr. 

For t < tr, < and, hence, increases. The evolution can then 
be analysed along similar lines as above. In particular, if x remains in the 
range < x < 1 then no singularities arise. However, depending on the 
details of matter content and the values of e"^ and x niay oscillate back 

and forth. But, the singularities can potentially arise only when x ^ 0, or 
X — > 1 where f2 — > oo. x — > is analysed above and x 1 will be analysed 
presently. 

3c. t >ti : A{ti) and xiU) are positive, and the evolution will now lead 
to the present day universe. For t > ti, GUT symmetry breaking, inflation, 
and other (matter) model dependent phenomena may occur, but are all non 
singular. Also, in all such phenomena, both and x increase. Thus, as t 
increases beyond tj, x increases and eventually x 1- Consider now this 
limit, where singularities may arise because — > oo. 

Asx~^l) f^— — x)~^"; see (H). Then, to an excellent approxima- 
tion, cr(t) is constant and equation (0) can be solved exactly. Proceeding as 
in 3b above yields the unique solution, which also describes the present day 
universe: 

^A-Ao ^ ^ X = 1 - Xot~ , (11) 

where Aq and Xo > are constants. It can now be easily seen that the 
quantities in ( P^ are all finite, implying that all the curvature invariants are 
finite. Thus, there is no singularity as x 1- This shows that there is no 
singularity for t > U (and also for t < tr in 3b above). 

Thus, in the present day universe in our model, x 1 and i7 — > oo. Also, 
as follows from (||), ^^{today) oc (1 — x)^"^^ ~^ 0- Therefore, our model 
satisfies the observational constraints imposed by solar system experiments, 
viz. VLitoday) > 2000 and ^^{today) < 0.0002 (see g, pg. 117, 124-5, and 
339). ^ 

4. We have thus shown that a homogeneous isotropic universe, such as 
our observed one, evolves in the present model with no big bang or any other 
singularity. The time continues indefinitely into the past and the future, 
without encountering any singularity. This is a generic result valid for any 
function Q satisfying only the constraints (0). Note that the evolutions given 
in 1^, ^ for specific choices of fi(x) and/or p and 7 all conform to the evolution 
described here, and are singularity free only when fi(x) satisfies (0). 
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Now, a large class of fl satisfying (Q) can be easily constructed. An 
important question to ask, however, is whether such a function Q can arise 
from string theory. We now discuss critically a few possibilities. 

In this context, Einstein metric formulation given in is more natural 
where is the dilaton and x(0) the arbitrary function. The strong coupling 
limit in string theory, relevant for singularities, corresponds to diverging ef- 
fective Newton's constant (= j^^) and, thus, to the limit 0- this 

limit, — oo and Inx — > in our model (see equation (j^) and 0). 

(1) Note that, in string theory, ln(x(0)) can be thought of as Kahler 
potential for 0, which is expected to be modified at strong coupling by non 
perturbative effects. But, such modifications are exponential in nature and, 
hence, are unlikely to lead to a x(0) as required here. 

(2) However, at strong coupling, here may not be the stringy dilaton, 
but instead a combination of the stringy dilaton and other compactification 
dependent moduli. In that case, after writing the relevant effective action as 
given in |]^, the Kahler potential In^ for this 'effective dilaton' may well 
turn out to be of the required form. 

(3) Another, perhaps more promising, possibility is the following which 
arises when JIq = | (see equation (|^) and |^). This case corresponds to a five 
dimensional space time compactified to four dimensions on a circle [0. In 
recent developments in string theory at strong coupling, analogous phenom- 
ena relating a (i- dimensional theory and a (d + l)-dimensional theory on a 
circle are found to occur: Using S-duality symmetries, Witten has discovered 
that the ten dimensional string at strong coupling is an eleven dimensional 
(M-)theory compactified on a circle [0, A similar phenomenon, that the 
three dimensional string at strong coupling is a four dimensional theory com- 
pactified on a circle, is at the heart of Witten's novel proposal for solving 
the cosmological constant problem |T^; for its possible stringy realisation see 
[Q. A similar phenomenon in four /five dimensions, if exists, is likely to lead 
to a x(0) as required here. 

Admittedly, at present, these are plausibility arguments only. Neverthe- 
less, given the elegant way the present model resolves the big bang singularity 
it is worthwhile to derive it, perhaps along the above lines, from a funda- 
mental theory such as string theory. On the other hand, however, such a 
model, even if not derivable from string theory, can stand on its own as an 
interesting model for singularity free cosmology. One may then study its im- 
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plications in other cosmological and astrophysical contexts, which are hkely 
to be novel |0. Such studies are fruitful and are likely to provide valuable 
insights. 
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Appendix: Curvature invariants are constructed using metric tensor, 
Riemann tensor, and covariant derivatives. When the metric is diagonal, 
every term in any curvature invariant can be grouped into three types of 
factors: (A) J^^^^g^^^^g^ R^.xr. (B) V¥^W^ ^"^^ (^) ^ 



acting multiply on (A) and (B) type factors (no summation over repeated 
indices) . 

Using the metric ds^ = — dt^ + e^^^*^ i i-^L'^ +r^(i5'|) and evaluating explic- 
itly the above factors, one obtains that (A) and (B) type factors are functions 
of A, A, and e~^. The action of (C) produces extra time derivatives. Thus, 
it follows that the curvature invariants are functions of and n > 1. 

By repeated use of equations (§)-(§), for any n > 1 and, hence, any 
curvature invariant, can be expressed in terms of the following quantities: 

The resultant expressions are finite if the above quantities are finite. The 
algebra is straightforward and, hence, we omit the details here. 

Therefore, a sufficient condition for all the curvature invariants to be 




finite is that the quantities in ^TW) be all finite. 
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